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The interaction Hamiltonian for a system of polarons a la 
Feynman in the presence of long range Coulomb interaction 
is derived and the dielectric function is computed in mean 
field. For large enough concentration a liquid of such parti- 
cles becomes unstable. The onset of the instability is signaled 
by the softening of a collective optical mode in which all elec- 
trons oscillate in phase in their respective self-trapping poten- 
tial. We associate the instability with a metallization of the 
system. Optical experiments in slightly doped cuprates and 
doped nickelates are analyzed within this theory. We discuss 
why doped cuprates matallize whereas nickelates do not. 
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An electron moving in a highly polarizable lattice can 
distort the environment and create a potential which 
self-traps the electron. The resulting object can move 
by dragging the distortion resulting in a heavy particle 
called polaron Q|. Apart from its motion as a whole 
the polaron has an extra degree of freedom which can 
be thought as the electron oscillating in the self-trapping 
potential. Because of this internal mode a polaron can 
be polarized in an external electric field. 

In this Letter a model Hamiltonian that takes into ac- 
count the interactions due to mutual polarization is de- 
rived. It is shown that a polaron liquid or more generally 
a liquid of polarizable particles becomes unstable for high 
enough concentrations due to the dipole-dipole interac- 
tion. The onset of the instability is indicated by a soft- 
ening of the internal mode of the polaron which can 
be measured in optical experiments |^-^ . Recent exper- 
imental results in cuprates and nickelates are analyzed 
within this scenario (^-|^. 

For simplicity the present treatment is semiclassical 
(in the spirit of the Drude model ^) however we expect 
that the same effect will show up in a fully quantum 
mechanical computation. 

We concentrate on the case of dielectric polarons but 
we believe our results are valid also for other kinds of 
heavily dressed particles in solids as long as they have a 
polarizable internal degree of freedom. 

Consider an electron in a dielectric with static dielec- 
tric constant eg dominated by lattice polarization. If the 
electron is not allowed to move it will generate a radial 
polarization of the lattice towards itself. This distortion 
implies that some ionic positive charge (q) will accumu- 



late in the vicinity of the electron. The charge induced 
in the lattice can be estimated as 



q 
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Now let us free the electron. If the electron binds to 
the induced lattice charge we have a dielectric polaron. 
In this case for the electron to move it has to drag the 
distortion of the lattice. 

To describe this situation we use Feynman's effective 
model for the polaron The effect of the distor- 

tion is mimic by a heavy fictitious particle of mass M. 
The electron, of mass to, is coupled to the heavy parti- 
cle with a harmonic potential. The result is a composite 
particle. The non interacting Hamiltonian for a liquid of 
such particles reads Hq = ^ • Hi with 
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The first term in Hi is the kinetic energy of the electron. 
The second term is the kinetic energy of the heavy mass 
representing the surrounding deformation. The third 
term is the coupling between the electron and the dis- 
tortion. Pi, Ti {Pi, Ri) are the momentum and position 
of the electron (fictitious mass) i. 

In the original formulation by Feynman the action cor- 
responding to Eq. (^) is regarded as an effective action for 
the dynamics of a polaron. The parameters are obtained 
variationally from the Frohlich model of a dielectric 
polaron. Hare we take as granted that Hq is a first ap- 
proximation to describe the dynamics of a polaron liquid 
and explore the consequences of adding a Coulomb term. 
On physical grounds we will assume the charge q is dis- 
tributed with spherical symmetry centered on Ri. 

To derive the interaction Hamiltonian consider two 
fixed electrons in the dielectric at ri,r2. As before each 
electron will induce a close-by charge in the lattice of 
magnitude +q. The total electrostatic energy can be 
computed as usual. One gets self-energy terms plus the 
interaction term: 
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This energy includes the vacuum electrostatic energy be- 
tween charges plus the elastic energy stored in the dielec- 
tric: 
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The former is: 



Coul 



E, 



la- 



eSri -^2| 



(4) 



(5) 



i.e. the electrostatic energy between the screened 
charges. Eq. ffl) can be solved for the elastic energy: 
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The factors with the dielectric constant in Eq. (||) ensure 
that the elastic energy is zero for an infinitely rigid lattice 
(eo —>■ 1) or an infinitely soft lattice (eg — > oo). 

Suppose now that we move the electrons from the equi- 
librium position keeping ions fixed so that the distortion 
does not move. i.e. ^ Ri- Since the elastic energy de- 
pends only on the configuration of the lattice we should 
replace in Eq. (^) Ri. So that in general the elastic 
energy is: 
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where we eliminated eo using the definition of the induced 
lattice charge q [Eq. (^J. 

The Coulomb energy can be decomposed into ele- 
mentary Coulomb interactions. For general position of 
charges and distortions it reads: 
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which reduces to Eq. (|^) for = Ri. 

Adding again the elastic and Coulomb energy we can 
write the interaction Hamiltonian for a liquid of such 
particles: 
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where indexes i,j run over particles. 

A similar argument can be used to derive the inter- 
action Hamiltonian with an external electric field. The 
Coulomb contribution of the induced charge cancels with 
the elastic part and one obtains, 
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The many-particle Hamiltonian is H = Hq + H\^^ + He- 
It is convenient to change to center of mass variables, 
Pi = {rim + RiM)/{M + to), Ui = ri — Ri. Making 
a Taylor expansion for small Ui we obtain the following 
interacting Hamiltonian in the dipole approximation. 
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where we defined the matrix (f)^^{p) = 6^^/ p'^ ~3p^p,y/ 
and /i, ly are Cartesian indexes. 

Hamiltonian Eq. (||) have the essential ingredients to 
treat the interacting polaron problem at concentrations 
such that the interparticle distance exceeds the polaron 
radius. 

Now we compute the dielectric constant of such a sys- 
tem. To obtain equations of motion for one particle 
in the mean field of the others we compute the forces 
= —dH/dui and Fp. = —dH/dpi. The former is 
the force that polarizes particles and the latter is the force 
acting on the center of mass. The equations of motion in 
mean field read, 

nil = —ku + AnLnequ — tttE (12) 
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where u =< Ui >, etc. n is the density of particles, 
l//i = 1/to-|-1/M, and L is a geometric factor discussed 
below. 

A crucial point in the derivation is the evaluation of 
the force due to dipole-dipole interactions Fl = equ. < 
'^jj^i'PiPi ~ Pj) >• This is the well known Lorentz- 
Lorenz local field i.e. the dipolar field at pi due to dipoles 
at positions pj averaged over the position of the dipoles. 
Using well known results from the theory of dielectrics 
[|l| we obtain the second term in Eq. (^2|). Either for a 
random isotropic distribution of dipoles or for a cubic ar- 
ray L = 1/3. For other distributions L can be computed 
using the results of Ref . ||l^ (see below) . 

By putting E = E^e^'^'^* we solve for the dipole mo- 
ment n and compute the polarization vector P = nil 
where in the original variables the dipole moment is 
Hi = -eri+ qRi . 

Finally from the relation between P and E one ob- 
tains the following dielectric function for the interacting 
polaron liquid 
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where we have introduced phenomenological inverse re- 
laxation times 7, 7' and defined the bare plasma fre- 
quency ll>p = Aime^ /m and the polaron plasma frequency 
Vil = 47rne(e - q)/{M + m) = 47rne^/eo(Af -|- m). In the 
latter M + m can be identified with the polaron effective 
mass. We have also defined the renormalized frequency 
of the internal mode 
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where ujq = \/k/ii is the frequency of the internal mode 
of an isolated polaron. 

The dielectric function Eq. (|l^) takes a very simple 
form: the sum of a Drude term plus a collective ex- 
citation at frequency Wcoii- The latter consists of all 
electrons oscillating in phase inside their respective self- 
trapping potentials. In this oscillation polarons develop 
time-dependent in-phase dipole moments. The dipole- 
dipole interaction tends to soften this mode and the effect 
increases as interparticle distance decreases i.e. density 
increases (oc uip). Eq (|lj) can be seen as a generaliza- 
tion of both Drude model and Clausius-Mossotti equation 
0. Notice that the f-sum rule is exactly satisfied. 

At some critical density given by ATrric/e^m = 
ujQe/qL{l + m/M) the energy cost to displace the elec- 
trons away from the respective lattice distortion vanishes 
and an instability occurs. Within this simplified model 
we cannot describe the new phase that arises for n > Uc- 
In fact we assumed above that electron and distortion 
are bound and, even more, that the binding potential is 
harmonic. Both approximations will break down close to 
Uc. Interestingly, if a positive quartic term were present 
in the distortion-electron potential, the system would be- 
come a liquid crystal of ferroelectric polarons above the 
critical density. 

A less exotic possibility is that polarons start to dis- 
sociate. This can occur either abruptly or in a continu- 
ous way. In the former case all polarons collapse at Uc 
whereas the latter case can be realized by a two com- 
ponent system of coexisting polarons and free electrons 

Recently the softening of a polaron band has been ob- 
served as a function of doping 1^1 in the Nd2Cu04_j, sys- 
tem. We have used Eqs. (|l^),(|l5|) to analyze those data. 

Slightly doped cuprates show a polaronic band, phonon 
bands, a mid-IR band and a charge transfer band 
(we neglect much weaker magnetic bands [^,^). Fol- 
lowing Calvani and collaborators we fitted reflectivity 
data (|,^,|| on Nd2Cu04_j^ with a Drude-Lorentz model 
for the dielectric constant (a sum of Lorentzians plus a 
Drude term, see Ref. The model dielectric function 



is of the same form as in Eq. (14) but for the addition 
of phonons and of high energy electronic contributions, 
and with an appropriate eoo replacing 1. Because of the 
electronic screening ujp and flp in Eq. (|l^) should be con- 
sidered screened plasma frequencies. Two Lorentzians 
were used above 10^ cm~^ to fit the charge transfer band 
and higher energy contributions. Four Lorentzians were 
used below 800 cm^^ to model the TO phonons, one 
Lorentzian above 5000 cm^^ (depending on sample) 
were used to fit the mid-IR band (MIR) and finally one 
Lorentzian below 2000 cm^^ were used to fit the inter- 
nal mode of polarons. In Fig. (|l|) we show the result- 
ing optical conductivity 17(0;) = lme{uj)uj / Air , excluding 



the electronic contributions and the phonons. Neither n 
nor y are well controlled variables so we use Wp(oc n) as 
our control parameter which can be determined directly 
by adding the polaron and the Drude oscillator strength 



[Eq. (14)]. In the inset of Fig. |^ we show uj^^n vs. ujp. 



Both quantities are obtained from the fits. 

As found by Lupi et al. ^ in this doping range we see a 
rapid decrease of the internal mode energy. To estimate 
the rate of decrease we neglect the higher doping point (it 
is in the region where our approximations break down) 
and we do a linear regression. We obtain a slope of —0.5. 
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FIG. 1. Optical conductivity obtained from the fit of the 
model dielectric function to reflectivity data at 300K 

as explained in the text. Only Drude and internal mode po- 
laron contributions have been included. The upper curves 
have larger doping y or equivalently larger density of parti- 
cles n. The area under the curves is proportional to uj^. The 
maximum of the broad structure approximately corresponds 
to Wcoii. Inset: The parameter w^oii vs. uj^ from the fits. The 
error bars refiect uncertainties in the fit. The line is a linear 
regression excluding the higher doping point. 



To get a theoretical estimate of the slope in Eq. (15) we 
need the Lorentz-Lorenz local field factor L. Using the 
results of Ref. |jl^ we estimate L — .6 for Nd2Cu04_y 
and for the electric field of the radiation parallel to the 
Cu-0 plane. 

Also the induced charge in the lattice will not be given 
any more by Eq. (|^). In fact Eq. (|l|) gives the total in- 
duced charge but part of this charge which we call 
is induced in the electronic degrees of freedom. This 
should be subtracted from +q and added to — e since 
it will follow antiadiabatically the electron and hence is 
positioned at r^, not at Ri. This effect is taken into ac- 
count by replacing q/e hy {q — q^) / {e — q^) in Eq. (|T^). 
qe can be estimated, using long wave length Lyddane- 
Sachs-Teller like arguments [|l|, to be of the order of 
qe = e(l - l/eoo)(l + 2/eo)/(l + Using eq = 30 

and = 5 [H one gets {q - qe)/{e - qe) = 0.91. Finally, 
as a first approximation, we can take (1 -I- M/m) ~ 1 
so we obtain for the slope -0.55 in excellent agreement 
with the experiment. The linear regression extrapolates 
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to Wcoii = at a critical value w^^^ = 7.6 x 10^ cm"^. 

For n > Tie a. sudden increase of Drude weight is ob- 
served (not shown), that we associate with metallization. 
This occur at much lower doping than the insulator su- 
perconducting transition observed at lower temperatures 
pO| . We believe that more complicated phenomena like 
charge ordering may change the picture at low temper- 
atures |2|. Quantum corrections will also be important 
close to Uc and at low temperatures. 

Lupi et al. find that the collective mode subsists at a 
small finite frequency for n > ric and then monotonously 
decreases as doping is increased at a much slower rate . 
This suggests a scenario where the polarons stabilize at a 
concentration slightly below ric and added carriers go to 
free states. Different versions of such two-fluid model has 
been considered in the literature |14-I^]. Interestingly, 
this coexistence implies that for a finite doping range the 
system is at the verge of a dielectric anomaly with a soft 
electronic mode. At larger doping, polarons are expected 
to become unstable due to short-range interactions and 
screening and a first order transition will occur to a one 
component Fermi liquid. We speculate that this scenario 
can explain the phase diagram of cuprates. In fact the 
soft collective mode can explain many of the anomalous 
Fermi liquid properties pT[ | that these systems show at 
optimum doping levels including pairing formation and 
superconductivity due to exchange of the soft mode bo- 
son. At overdoping the experimental optical conductivity 
shows a single component Drude behavior which suggest 
a first order transition |^^. A related scenario includ- 
ing the softening of an optical electronic mode associated 
with matallization has been proposed previously . 

The present theory also explains why a material like 
doped La2Ni04 never becomes metallic. In this system 
the observed energy of the internal mode for small doping 
is roughly by a factor of 3 larger than in Nd2Cu04_j^ . 
Then if we assume that the slopes of w^^jj vs. uu'i are sim- 



ilar in the two systems we expect, according to Eq. (15), 
that the critical plasma frequency to^ ^ in the nickelate is 
by a factor of 9 larger than in the cuprate. On the other 
hand, the observed spectral weight in this system reaches 
a maximum and decreases again as a function of doping 
We estimate the maximum by integrating the 
fitted line shapes in Ref. to be roughly 9 x 10^cm~^ 
which is almost by an order of magnitude smaller than 
the estimated ^. This means that on undimensional- 

ized scales (w^ou/'^o ^p/^px) ^'^^y ^'^^^ ^ -'^^■^ 
of the curve in the inset of Fig. H is physically accessible 
and the system never reaches the point in which it should 
metallize. Meaningfully a modest softening of the inter- 
nal mode is actually observed [^,^ in this relatively small 
^p/^p.c variation range consistent with Eq. (|l5|). 

To conclude we have derived a many body Hamiltonian 
that takes into account dipole-dipole interactions in a po- 
laron system in a general way. We obtained the dielec- 
tric constant in mean field for a liquid of such particles 



and showed that the system becomes unstable for large 
concentration. We argue that the theory explains a rapid 
softening of the internal mode of the polaron recently ob- 
served in cuprates ||^ and discussed the doping-induced 
matallization of cuprates and nickelates. We speculate 
that the same framework can be used to explain the oc- 
currence of matallization or not as a function of doping 
in a wide range of materials where the electron phonon 
interaction dominates. 

We thanks P. Calvani, S. Lupi, M. Capizzi, S. Fratini 
and P. Quemerais for enlightening discussions and Lupi 
et al. for providing us their reflectivity data before pub- 
lication. 
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